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P.  R.  Fischer,  J.  L,  Daschbach  and  G.L.  Richmond 
Dept,  of  Chemistry 
University  of  Oregon 
Eugene,  OR  97403 


Abstract: 

The  optical  second  harmonic  (SH)  response  from  Si(lll)  electrode  surfaces  has 
been  studied  and  has  been  found  to  be  highly  potential  dependent.  For  both  H-terminated 
Si(lll)  surfaces  in  NH4F,  and  oxide  covered  surfaces  biased  in  H2SO4,  the  phase  and  the 
amplitude  display  a  potential  dependence  which  we  attribute  to  field  effects  within  the 
semiconductor  space-charge  region  and  at  the  surface  of  the  Si(l  1 1)  electrode.  The  studies 
arc  the  first  to  demonstrate  that  for  Si(l  1  l)/oxide  samples,  the  potential  dependence  in  the 
SH  phase  can  be  correlated  with  oxide  thickness,  as  demonstrated  by  examining  samples  of 
0-40  A  thickness. 
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I.  INTRODUCTION 
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Although  Si/electrolyte  and  Si/Si02/electrolyte  interfaces  have  been  extensively 
investigated,  many  issues  regarding  the  structural,  electronic  and  electrochemical 
properties  of  these  interfaces  remain  of  interest.  Understanding  the  nature  of  these 
interfaces  has  much  technological  relevance.  In  recent  years  significant  progress  has  been 
made  in  characterizing  the  Si/electrolyte  interface  with  FT-IR  experiments,*’  ^  UHV 
transfer  experiments,^  and  STM^>  6,  7  an(j  AFM^  imaging  by  taking  advantage  of  the 
relative  stability  of  Si  surfaces  prepared  in  a  hydrogen  terminated  state.  However,  these 
probes  are  inapplicable  for  measuring  the  Si/Si02  interface  once  several  overlayers  of  Si02 
are  present.  In  this  letter  we  report  the  results  of  the  first  surface  second  harmonic 
generation  (SHG)  study  which  explicitly  probes  the  properties  of  the  Si/SiC>2  electrolyte 
interface  as  a  function  of  oxide  thickness.  The  results  are  compared  with  the  potential 
dependent  response  for  H-terminated  Si  surfaces.  Because  of  the  inherent  surface 
sensitivity  of  SHG,  it  is  a  useful  probe  of  interfacial  regions,  particularly  buried  interfaces 
which  are  difficult  to  probe  by  other  techniques.^  Although  SHG  and  sum  frequency 
generation  (SFG)  have  been  used  in  many  studies  to  examine  Si(l  1 1)  in  UHV*0> 1  *  and  in 
a^,  12,  13  application  to  Si/electrolyte  and  Si/SiC>2/electrolyte  interfaces  have  been  almost 
nonexistent.  Of  the  two  earlier  electrochemical  studies  reported  in  this  area,*^>  15  both 
were  performed  prior  to  recent  advances  in  surface  preparation  procedures  and  neither 
have  involved  the  polarization  and  surface  orientation  experiments  necessary  to  understand 
the  origin  of  the  potential  dependent  response  as  has  been  done  here. 

The  most  striking  result  is  the  strong  potential  dependence  that  we  observe  in  the 
relative  phase  of  the  SH  response,  as  measured  with  rotational  anisotropy  experiments,^ 
from  both  the  H-terminated  surfaces  and  surfaces  with  oxides.  We  attribute  the  potential 
dependence  for  both  the  predominately  H-terminated  and  the  oxide  covered  surface  to  the 
change  in  the  static  field  at  the  surface  and  within  the  space-charge  region  (SCR)  of  the 
semiconductor.  In  the  presence  of  the  insulating  oxide  layer,  the  potential  dependence  can 
be  interpreted  in  terms  of  a  screening  of  the  portion  of  the  applied  dc  field.  The  work 
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demonstrates  that  such  phase  measurements  provide  a  new  means  of  characterizing  the 
properties  of  the  space  charge  region  and  the  field  drop  across  the  insulating  oxide  layer. 
The  ability  to  observe  potential  dependence  and  its  origin  has  been  questionable  up  to  this 
point  because  of  the  relatively  large  bulk  response  from  semiconductors.  A  recent  report 
has  investigated  the  effect  of  the  depletion  layer  electric  field  on  the  SHG  response  from 
GaAs(OOl)  under  UHV  conditions  but  this  study  was  limited  to  a  smaller  range  of  electric 
fields  than  reported  here  and  was  performed  only  on  native  surfaces.^ 

II.  THEORETICAL  CONSIDERATIONS 

When  a  potential  is  applied  across  a  semiconductor/electrolyte  interface,  the  second 
harmonic  response  can  be  written  as 

II  F  (&2,+  X?  +  S’A<D)  ||2  (1) 

where  F  corresponds  to  the  linear  Fresnel  factors,  X,T  and  X^'  are  the  surface  and  bulk 
susceptibilities,  is  the  effective  cubic  nonlinearity  arising  from  the  static  field 

which  includes  all  terms  that  vary  linearly  with  field  strength.  ^  A<X>  is  the  potential  drop 
across  the  semiconductor  space  charge  region  and  is  proportional  to  the  difference  between 
the  applied  (Eapp)  and  the  flatband  (Eft,)  potential.  When  the  SH  response  is  dominated  by 
the  cubic  nonlinearity  term,  the  observed  potential  dependence  should  be  parabolic  with  a 
minimum  near  the  flatband  potential.  Such  behavior  has  been  observed  in  numerous 
studies  of  metal/electrolyte  systems.  Under  conditions  where  the  surface  and/or  bulk 
quadratic  nonlinearities  dominate,  parabolic  potential  dependence  with  a  minimum  shifted 
from  flatband  would  be  observed. 

By  judicious  choice  of  input  and  output  polarizations,  and  by  monitoring  the 

variation  in  the  SH  response  as  the  cubic  sample  is  rotated  azimuthally,  a  detailed  analysis 
of  the  relative  contribution  from  the  various  tensor  elements,  X^  can  be  performed.  ^ 

The  observed  angular  variation  (rotational  anisotropy)  can  be  related  to  the  overall 
symmetry  of  the  surface  atomic  structure.  The  intensity  of  this  azimuthal  dependence  for  a 


(Ill)  surface  and  the  input  and  output  polarizations  used  here  can  be  described  by  the 
phenomenological  equations  ^ 

I(2ca)pp  «  I  a  (-)  +  c^cos(3<(>)  I2  (2) 

I(2<o)p,s  oc  I  b(3  ')sin(30)  I2  (3) 

where  the  azimuthal  angle  <j>  is  defined  as  the  angle  between  the  [211]  direction  and  the 
projection  of  the  incident  wavevector  parallel  to  the  surface.  The  subscripts  p  and  s  denote 
the  beam  polarizations  for  the  fundamental  and  second  harmonic  respectively.  The 
isotropic  coefficient  a  (°°),  and  the  anisotropic  coefficients,  b(3>  and  d3K  contain  both  the 
Fresnel  coefficients  as  well  as  the  material  susceptibility  elements.  The  isotropic 
coefficient,  often  discussed  as  the  out-of-plane  response,  contains  the  three  independent 
surface  tensor  elements  X*'  and  Xm  (i=x,y)  as  well  as  the  bulk  elements  C,  and  y. 
The  anisotropic  coefficients  contain  the  in-plane  Xm  surface  tensor  element  and  bulk 
element  Because  the  experiments  conducted  under  both  p-input  and  output  polarizations 
(Eq.  2)  involve  an  interference  between  a(°°)and  the  fits  to  the  data  provide  information 

about  the  relative  magnitude  and  phase  of  these  two  terms.  For  the  potential  dependent 
measurements,  X^y  has  the  same  form  as  xfjjj*  as  a  result  of  the  applied  field  being  parallel 

to  the  surface  normal  about  which  the  sample  is  rotated. 

IH,  EXPERIMENTAL  METHODS 

The  optical  measurements  employ  the  fundamental  output  from  a  10  Hz  Nd:YAG 
laser  producing  10  nsec  pulses  which  illuminate  the  surface  at  a  32°  incident  angle.  Pulses 
of  =0.3  J/cm2  were  used  which  is  below  the  damage  threshold  for  silicon.  The  optics, 
detection  instrumentation  and  electrochemical  cell  are  similar  to  that  described 
previously.^®  The  silicon  wafers  are  n-doped  with  phosphorus  with  a  resistivity  of  3.0  to 
6.5  Q-cm  which  represents  a  defect  doping  density  of  =1015  cm-3.  The  0.66-0.71  mm 
thick  samples  are  degreased  by  ultrasonification  in  separate  baths  of  methylene  chloride, 
acetone,  and  methanol  and  dried  with  nitrogen.  The  back  of  the  wafer  is  etched  for  one 
minute  in  48%  hydrofluoric  acid  to  remove  the  native  oxide  and  then  mounted  on  Ga-In 


eutectic  that  has  been  placed  on  an  embedded  copper  contact  in  a  Kel-F  shaft.  A  mask 
containing  an  embedded  acid  resistant  fluorocarbon  o-ring  is  used  to  seal  the  surface  from 
the  electrolyte. 

The  silicon  wafer  is  prepared  by  etching  in  buffered  NH4F  solution  which  is  known 
to  leave  the  surface  in  a  H-terminated  stated  This  is  followed  by  immersion  in  0.1  M 
NH4F  if  a  hydrogen  terminated  surface  is  to  be  maintained,  or  in  0.1  M  H2SO4  if  an  oxide 
is  to  be  grown.  Oxides  are  grown  on  the  Si  samples  photoanodically  by  the  stepwise 
increase  in  potential  from  the  flatband  to  a  final  potential  of  +5.0  V  while  keeping  the 
anodic  current  below  30|iA/cm2.  A  HeNe  laser  is  used  for  illumination.  Thicknesses  of 
the  oxides  are  determined  by  a  combination  of  ellipsometry,  photocurrent  measurements, 
and  etch-back  times.  All  potentials  are  referenced  to  the  saturated  calomel  electrode 
(SCE).  Flatband  potentials  are  determined  by  photocurrent  transients. 

IV.  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  rotational  anisotropy  in  the  SH  response  with  p,p  polarization 
for  a  Si(lll)  surface  studied  in  NH4F  (Fig.  1(a))  and  H2SO4  (Fig.  1(b))  at  potentials  near 
flatband,  -0.65  V  and  -0.6  V  respectively.  Both  samples  were  etched  for  3  min.  in  2.0  M 
NH4F  at  pH  4.5  to  produce  a  H-terminated  surface  prior  to  immersion.  The  three-fold 
symmetry  expected  of  this  (111)  surface  of  a  cubic  lattice  is  observed  in  both  solutions. 
The  results  for  the  surface  in  NH4F  gives  a  phase  angle  for  dtyd00)  of  25  ±  5°.  Previous 
studies  suggest  that  this  surface  in  NH4F  should  remain  predominately  H-terminated  when 
potentiostated  near  the  flatband.  It  should  also  be  relatively  free  of  any  photogenerated 
surface  oxides  since  in  fluoride  containing  electrolytes,  the  dissolution  of  this  oxide  should 
compete  with  the  photo-oxidation  process.  The  response  from  the  same  surface  studied  in 
H2SO4  is  similar  but  the  phase  angle  fluctuates  between  30°  -  60°.  The  flatband  potential 
also  varies  slightly  in  time.  We  attribute  both  variations  to  photogenerated  surface  species 
which  can  store  charge  at  the  interface.  From  photocurrent  measurements  we  estimate  the 
amount  of  photogenerated  species  to  be  less  than  3  ML. 


The  potential  dependence  of  a  Si(ll  1)  surface  in  0.1  M  NH4F  has  been  examined 
between  -0.65  V  (flatband)  and  +1.0  V.  The  experiments  have  been  conducted  at  <j)  =  30° 
so  as  to  isolate  the  potential  dependence  in  the  isotropic  and  anisotropic  contributions.  As 
shown  in  Figure  2(a),  both  the  isotropic  and  anisotropic  response  are  found  to  have  a 
parabolic  potential  dependence  with  a  higher  overall  signal  level  from  the  latter.  The 
isotropic  response  which  would  most  readily  couple  to  the  applied  static  field  has  a 
minimum  near  +0.26  V,  =900  mV  from  the  flatband  potential.  This  field  can  couple  to  a 
depth  of  the  space  charge  region,  which  for  this  surface  biased  at  +0.3  V  is  on  the  order  of 
1200  nm.  The  anisotropic  in-plane  response  has  a  minimum  near  +1.6V,  even  further 
away  from  flatband.  In  H2SO4  the  minima  for  the  isotropic  and  anisotropic  responses  are 
also  shifted  from  flatband. 

The  fact  that  the  minimum  signal  does  not  occur  at  the  flatband  for  the  isotropic  or 
anisotropic  response  in  either  solution  suggests  that  terms  other  than  x?J  Ad>  in  Eq.  1  must 

be  contributing.  xT  *s  the  likely  factor  considering  that  the  penetration  depth  of  the  light 
is  on  the  order  of  I  cm,  considerably  deeper  than  the  SCR.  Previous  studies  in  air  have 
determined  that  the  SH  contribution  from  the  surface  and  the  bulk  of  Si(l  1 1)  are  of  similar 
magnitude.  *3  Further  evidence  for  the  importance  of  the  bulk  response  comes  from  the 
observation  that  when  Si(lll)  surfaces  are  roughened  by  etching  in  48%  HF  prior  to 
introduction  into  the  NH4F,  the  rotational  anisotropies  for  the  roughened  and  unroughened 
samples  are  quite  similar.  'iTiere  is  no  evidence  of  an  isotropic  response  from  a  disordered 

silicon  surface  superimposed  on  the  response  from  the  crystalline  lattice.  The  intensity  of 
the  overall  response  is  however  slightly  larger  in  the  roughened  case,  suggesting  that  jj^2) 

is  not  insignificant  Analysis  of  the  response  in  Fig.  2(a)  based  on  Eq.  1  indicates  that  the 
sum  of  (*<2)  +  X?)  is  at  least  of  comparable  magnitude  to  the  contribution  from 

Si**  • 

In  earlier  work  it  was  assumed  that  the  electrode  polarization  was  absolute  with 
respect  to  zero  applied  potential  (flatband)  and  a  model  of  quadratic  dependence  of  the 
SHG  with  applied  potential  was  found  to  fit  the  data.  ^  ^  Although  the  flatband  was  not 
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measured  in  these  earlier  studies,  the  current- voltage  data  presented  suggests  that  the 
flatband  potential  was  well  negative  of  0  V  (vs.  SCE).  With  this  correction  the  data  clearly 
do  not  fit  a  quadratic  dependence  in  applied  field  strength  but  does  resemble  the  data 
reported  here. 

In  a  similar  manner  the  SH  response  has  been  studied  for  Si(lll)  surfaces  for 
which  varying  thicknesses  of  Si02  have  been  photoanodically  grown  on  the  surface. 
Figure  3  shows  the  rotational  anisotropy  from  an  oxide  covered  surface  in  H2SO4  at  two 
potentials,  0.0  V  (Fig.  3(a)),  the  flatband  potential  at  this  oxide  thickness,  and  +3.0V  (Fig. 
3(b)).  The  potential  "window"  for  study  of  the  oxide  coated  surfaces  is  wider  due  to  the 

o 

insulating  nature  of  the  overlayer.  The  oxide  layer  is  estimated  to  be  =25  A  based  on  etch 
back  times ^  and  ellipsometry  measurements.  Because  the  electrochemical ly  grown 
oxides  have  substantial  water  content  in  the  first  monolayers  of  oxide  with  this  ratio 
decreasing  as  the  thickness  increases,  we  view  our  ellipsometrically  determined 
thicknesses  as  an  upper  limit  In  comparing  the  response  from  the  Si(l  1 1)  surfaces  of  Fig. 
1(b)  measured  at  flatband  with  the  oxidized  sample  in  Fig.  3(a),  two  important 
observations  can  be  made.  The  rotational  anisotropy  persists  in  the  presence  of  the  oxide 
and  the  signal  from  the  oxide  covered  surface  is  enhanced  relative  to  the  more  oxide  free 
sample  in  Fig.  1(b).  The  former  observation  is  not  surprising  since  the  SH  response  from 
the  amorphous  oxide  overlayer  is  very  small, The  signal  enhancement  observed  in  the 
presence  of  oxide  demonstrates  that  the  SHG  is  sensitive  to  the  Si(ll  1)  surface  adjacent  to 
the  oxide  and  is  not  merely  a  bulk  response.  This  enhanced  signal  could  be  due  to  a 
crystalline  interfacial  region  between  the  Si  and  the  oxide  overlayer  as  suggested  in 
previous  studies  on  thermal  oxides.^  It  is  more  likely  that  the  effect  can  be  attributed  to 
an  interfacial  strain  in  the  Si  lattice  as  has  been  observed  in  a  strong  resonant  enhancement 
at  SH  photon  energy  of  3.3  eV  and  attributed  to  a  direct  transition  taking  place  in  a 
perpendicularly  strained  silicon  layer  at  the  thermal  oxide  Si/Si02  interface.^4 
Comparison  of  thermal  and  anodic  oxides  will  be  reported  elsewhere. 
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As  with  the  Si(l  1 1)  surfaces,  there  is  a  strong  potential  dependence  in  the  response 
from  the  Si(l  1  l)/SiC>2  interface  that  is  manifested  in  both  the  relative  magnitude  and  phase 
of  the  rotational  anisotropy  (Compare  Fig.  3(a)  and  3(b)).  For  an  applied  field  of  +3.0  V 
beyond  the  flatband  potential,  which  for  this  surface  is  0.0  V,  the  relative  phase  angle 
between  the  anisotropic  and  isotropic  response  changes  from  48°  to  129°.  The  relative 
magnitude  changes  by  a  factor  of  two.  The  measured  minima  in  the  anisotropic  and 
isotropic  response  occurs  near  +3.5V  and  +2.5  V  respectively,  both  minima  being  far  from 
the  flatband  potential.  As  with  the  Si(lll)  surfaces,  the  potential  dependence  for  the 
oxidized  surfaces  does  not  follow  the  predicted  form  if  Ad>  were  the  dominant  factor. 

To  understand  the  strong  potential  dependence  in  the  relative  phase  angle  observed 
in  Fig.  3  we  have  conducted  rotational  anisotropy  experiments  at  a  series  of  potentials  for 
surfaces  of  different  oxide  thicknesses.  A  portion  of  the  results  are  summarized  in  Table  1. 
Because  the  flatband  potential  varies  with  the  oxide  thickness,  all  data  is  referenced  to  the 
flatband  potential  of  each  sample  as  determined  by  photocurrent  transient  measurements. 
For  all  samples,  the  potential  dependence  follows  the  trend  observed  in  Figure  3  for  the 
25A  sample  where  the  phase  angle  increases  with  increased  positive  potential.  At 

flatband,  the  fits  to  the  anisotropies  for  the  different  oxide  covered  samples  show  little 
variation  as  one  would  expect  since  +  X^)  should  be  similar  for  samples  with 

more  than  a  few  monolayers  of  SiC>2.  The  fact  that  the  surface  in  NH4F  has  a  slightly 

different  magnitude  and  phase  (Fig.  1(a))  than  the  H2SO4  case  is  consistent  with  a  different 
X^  due  to  the  presence  of  oxidized  surfac ;  species.  At  potentials  positive  to  flatband,  a 

progressively  larger  Vapp  is  necessary  to  achieve  the  magnitude  of  the  phase  angle  that  is 
obtained  for  surfaces  with  thinner  oxide  coatings.  For  example,  the  phase  angles  obtained 

o 

for  the  40  A  sample  are  shifted  by  approximately  +1  V  relative  to  the  phase  angles 
obtained  from  the  25  A  sample.  We  attribute  the  potential  dependence  in  the  phase  angle 
to  the  variation  of  the  field  at  the  semiconductor  surface  and  within  the  space  charge 
region.  Because  of  the  insulating  nature  of  the  oxide,  the  field  at  the  Si(lll)/Si02 
interface  is  no  longer  Ad>  but  is  reduced  by  the  potential  drop  across  the  oxide  layer,  AV0x- 
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o 

This  would  suggest,  for  example,  that  the  additional  15  A  of  oxide  on  the  surface  screens 

o 

approximately  IV  of  the  overall  applied  field  relative  to  the  25  A  sample.  Later 
publications  will  address  in  more  detail  the  correlation  between  phase  and  amplitude 
changes  in  the  SH  response  with  oxide  thickness,  doping  density,  and  oxide  composition. 

In  summary,  SHG  studies  of  hydrogen  terminated  Si(lll)  electrode  surfaces  and 
Si(lll)/SiC>2  electrode  surfaces  of  varied  oxide  thickness  have  been  examined  to 
determine  the  origin  of  the  nonlinear  response.  The  potential  dependence  from  the  initially 
prepared  H-terminated  surfaces  in  NH4F  suggests  that  the  potential  independent  bulk 
response  is  significant  in  contrast  to  previous  studies^*  ^  which  have  described  the 
response  solely  in  terms  of  the  field  dependent  third  order  term.  For  oxidized  surfaces  the 
response  from  the  Si(lll)/SiC>2  interfacial  region  is  highly  sensitive  to  the  band  bending 
conditions  in  the  semiconductor  and  the  insulating  nature  of  the  oxide  overlayer.  The 
studies  show  the  sensitivity  of  SHG  to  the  electrostatic  properties  of  the  Si/SiC>2  electrolyte 
interface  and  demonstrate  the  utility  of  extending  this  method  to  semiconductor  electrode 
surfaces.  Studies  are  underway  to  provide  a  more  quantitative  comparison  between  oxide 
properties  and  the  potential  dependent  response  for  these  surfaces. 
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Figure  1.  Rotational  anistropy  from  Si(lll)  for  p-input  and  p-output  polarizations,  (a) 
Si(l  1 1)  in  0.1  M  NH4F  biased  at  the  flatband  potential  of  -0.65  V.  Best  fit  of  Eq.  3 
to  the  data  yields  =  2.2c*21*;  (b)  Si(lll)  in  0.1  M  I'2S04  biased  at  the 

flatband  potential  of  -0.6  V;  c^la°°  -  ITe*440. 

Figure  2.  Potential  dependence  in  the  isotropic  a(°°)  and  anisotropic  response  c<3>  for 
Si(l  1 1)  prepared  in  a  H-terminated  state  and  biased  in  0.1  M  NH4F.  The  isotropic 
response  was  monitored  with  p-polarized  incident  light  and  p-polarized  SHG  with  <|> 
=  30°.  The  anisotropic  response  was  monitored  with  p-polarized  incident  light  and 
s-polarized  output  light  at  <}>  =  30°. 

Figure  3.  Rotational  anisotropy  in  the  p-input  and  p-output  SH  response  from  an  oxidized 
Si(l  1 1)  sample  in  0.1  M  H2SO4  with  approximately  25  A  of  Si02-  (a)  At  0.0  V 
(flatband),  with  cMla™  =  2.3ei48‘;  (b)  At  +3.0  V,  c<3>Ja°°  =  4.7eil29\ 
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Table  1 :  Fits  to  Rotational  Anisotropy  Data:  dtyaf00) 


Oxide  Thickness 


V  app  *  V 

=  0  A  (H2SO4) 

15  A 

25  A 

40  A 

0  V  (FB) 

1.8  ei50* 

2.4  e*37’ 

2.3  e*48* 

2.5  e*42’ 

+0.5  V 

2.9  e155’ 

2.8  e*»* 

2.8  e^42* 

+1.0  V 

4.2  e'75’ 

3.4  ei72* 

3.2  e>51’ 

+1.5  V 

4.6  e'1^* 

4.0  ei89* 

3.5  e*8* 

+2.0  V 

4.0  eil26‘ 

5.0eiH3- 

3.7  e'70’ 

+2.5  V 

4.7  e'129’ 

4.5  e'92 

+3.0  V 

5.1  e'110' 

+3.5  V 

4.8  e«20* 
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